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Abstract 

We have calculated the 1"*" and 3"*", T=0 states in ^"^N and ^'^B . In a neutron-proton RPA 
model these two nuclei are described by the same set of equations. We first show that a bare 
Minnesota interaction leads to too weak binding in both nuclei. Furthermore it does not produce a 
3^ ground state in ^^B as it should. Including medium effects as an exchange of phonons between 
the neutron-proton pair cures the desagreement in ^^N but still gives a 1^ ground state in ^'^B with 
the 3"*" as an excited state. The same study with a Gogny effective interaction reproduces nicely 
the properties of both nuclei: same agreement in ^^N as previously when medium effeccts were 
introduced but now the 3+ in i°B becomes the ground state. This success suggests that through 
its density dependent term the Gogny interaction takes account of the presence of a three-body 
force which, in a shell model calculation, has been shown to be essential to give a 3^ ground state 
in . 
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I. INTRODUCTION. 



The T=0 neutron-proton pairing correlations in N=Z nuclei have motivated a number of 
theoretical and experimental works. In odd-odd N=Z nuclei, T=0 and T=l states are both 
present so that they offer a good place to make a comparative study of the two types of 
correlations. We have focused our attention on the two odd-odd nuclei of (N=Z=5) and 
-^^N (N=Z=7) for several reasons. First ^^N and ""^^B have a T=0 and a 3^ T=0 ground 
states respectively while the 0"*", T=l states are excited states and therefore suggesting 
stronger correlations in the T=0 channel. Furthermore because of charge invariance of the 
nucleon-nucleon interaction, the neutron-proton and neutron-neutron effective interactions 
should be the same in the T=l channel. Therefore an effective interaction and a nuclear 
model able to explain the O"*", T=l states in a core plus two neutrons system (like ^^C ) should 
also describe well the same T=l states in the core plus a neutron-proton pair (like ^^N ). 
We have then an opportunity to make a further test of our previous work in collaboration 
with J.C Pacheco on N=8 nuclei With the Gogny effective interaction or its zero 
range density dependent substitute jj], and in a two- neutron RPA model we were able to 
reproduce the 0+ states of ^^C -^°C, ^^Be-^Be and i^Li-^Li. The nucleus ^^N (^^B) described 
as a core of ^^C plus (minus) a neutron-proton pair is the analog of ^^C (^°C) described 
as a core of ^^C plus (minus) a neutron-neutron pair and our first aim is to check if the 
same interactions and the same RPA model where the two-neutron pair is replaced by a 
neutron-proton pair, will lead also to a good description of the 0"'' states in ^'^N and ^'^B . 

In the two-nucleon RPA model, successful in the description of the T=l states, we cal- 
culate the T=0, l"*" and 3"^ states in ^^N and ^°B . First we perform the calculation with a 
simple zero range force fitted to reproduce the deuteron binding energy. The results obtained 
with this simple potential show a too weak binding. Because a zero range force is a crude 
simplification of a realistic one, we have made the same calculation with the Minnesota bare 
interaction [S^ which has one short range component, two long range components and all 
exchange terms. The results are very close to the previous ones and show again a lack of 
binding. It is not surprising that bare interactions could not describe well bound many-body 
systems for which we know that medium effects are important. Medium effects on the bare 
nucleon- nucleon interaction may arise, at least partly, from the exchange of phonons between 
the two-nucleon pair. Indeed strong two-body correlations in nuclei manifest themselves as 
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very collective vibrationnal states at low excitation energies. These vibrations, or phonons, 
may be exchanged between the two nucleons and induce a modification of the nuclear inter- 
action. The presence of these phonons, besides their effect on the nuclear two-body force, 
has an other manifestation by modifying the average interaction of a single nucleon with the 
core due to the coupling between the nucleon and the phonons of the core. Such couplings 
modify the Hartree-Fock average potential and are so strong in nuclei like ^^Be and ^°Li that 
they are, at least partly, responsible for the inversion of the 1/2"^ and l/2~ states BBS. 
Both effects, on the two-nucleon force and on the nucleon-core interaction will be included 
in our two-particle RPA model for the description of ^^N and ^°B . 

The spectra obtained in ^^N and ^°B will be compared to those derived with the phe- 
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nomenological T=0 Gogny effective interaction [2|, |3| . From this comparison and comparison 
with shell model calculations using two-body and three-body interactions, we will try to 
understand better what is implicitly contained in this empirical force. 

In section II we present briefly the two-particle RPA model applied to a neutron-proton 
pair and precise our choice of single neutron and proton basis and our choice of two-body 
interactions. In section III we flrst present the results for the 0^, T=l states. Then we show 
the results obtained for the T=0 states with the two bare nucleon-nucleon interactions and 
discuss the contribution of phonon exchanges between the neutron-proton pair. Then we 
make the same study with the effective Gogny interaction. Comparing our results with those 
of shell model calculations we can give a qualitative interpretation of the medium effects 
contained in this very successful effective interaction. Section IV is devoted to our general 
conclusions. 

II. THE NEUTRON-PROTON RPA MODEL 

We describe ^^N and ^'^B as a core of ^^C in its ground state plus or minus a neutron- 
proton pair respectively and deflne two-body amplitudes as: 




(1) 
(2) 
(3) 




X, 
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Y- = {^'B\Aa\^''C) (4) 

where a,b. . . and a,/?,. . .represent configurations where the neutron and the proton are re- 
spectively in states unoccupied and occupied in the Hartree Fock ground state of the ^^C 
core. A^^^^ (Aa(a)) are operators which create (annihilate) a neutron-proton pair coupled to 
given spin and isospin (for simplicity we omit them in our equations). |^^C > represents the 
correlated ground state of the ^^C-core. Note that the Y-amphtudes, named the small RPA 
amplitudes, would be zero for an uncorrelated core. 

The two-nucleon RPA model has been described and used in a number of papers and we 
remind briefly that: 

-the same system of equations determines the amphtudes and energies of ^^N and 
which then are not independent. For a given spin and isospin these equations write as: 

(Q - ea)xa -Y.< > Xb-Y.< >xp^Q (5) 

b 13 

{n - ea)xa + X! < Oi\V\b > Xb + ^ < a\V\P > = (6) 

b 13 

where the eigenvalues,a;a and related to the amplitudes of eqs.(l-4) as explained 

below. The Ca and Ca are the unperturbed energies of the neutron-proton pair in states 
a and a respectively. The matrix elements of the neutron-proton interaction V have to be 
antisy mmetrised . 

-if the model subspace contains N configurations a,b,.. and M configurations a, (5 . . . , 
the RPA equations have N-l-M eigenstates with eigenvalues Q and eigenvectors x^ and x^. 
N of them correspond to ^"^N with: 



X 



(n) 



(7) 



and M to with: 



X-^'") = xt^ 



-a 



(8) 



y-(m) 



X 



(m) 



Eo(^^C ) is the ground state energy of the ^^C-core. The separation into two sets of solu- 
tions is unambigous and based on energy considerations and relative values of Xa and Xa 
amplitudes. 

- the amplitudes x^, and normalised according to: 



Ek'^r-EKr = 1 for ''N (9) 

a a 

= -1 for (10) 

-the only inputs of the calculation are the individual neutron and proton energies and 
the effective two-body interaction. 



A. Choice of the single particle basis 

We make a semi-phenomenogical approach to the RPA model. We replace the Hartree 
Fock average potential by a Saxon- Woods potential plus a spin-orbit force plus a phenomeno- 
logical surface potential fitted to reproduce the experimental single neutron energies in the 
field of ^^C and write the one neutron hamiltonian as: 

hn = tn + Vo (^/(r) - O.Urlils)^^^ + 6Vn (11) 

where: 

f(r) = [l + expC-^)]-' (12) 

with Vo=-50.5 MeV, a=0.75 fm, Ro=ro(12)i/3 with ro=1.27 fm. The last term SV^ is added 
to simulate medium effects due to the coupling of the neutron with the phonons of the core.. 
The shape of 6Vn is suggested by a semi-microscopic calculation of neutron- phonon couplings 
j?, [lO|. This contribution to the average one body potential depends on the neutron state 
and is written as: 

The coefficients a„ are fitted on the experimental neutron energies for lp3/2, lpi/2) 2s 
and ld5/2 states and put to zero for higher states. The proton energies are calculated by 
adding the Coulomb potential to the neutron hamiltonian of eq.(ll). In Table I are given 
the experimental energies and the corresponding calculated energies for the lowest, known, 
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lP3/2 

exp. cal. 


lPl/2 

exp. cal. 


2s 

exp. cal. 


1^5/2 

exp. cal. 


neutrons 
protons 


-18.72 -18.71 
-15.96 -15.54 


-4.95 -4.95 
-1.94 -1.98 


-1.85 -1.86 
0.42 0.25 


-1.1 -1.1 
1.61 1.35 



TABLE I: Experimental and calculated individual energies (in MeV) for neutrons and protons. 

neutrons and protons states. The results presented below have been obtained with the lp3/2 
proton energy replaced by the experimental value. However it has no effect on the ^^N 
spectrum and in gives both 1"*" and 3"*" energies are lowered by about 0.5 MeV what does 
not change qualitatively our discussion of results. The effect of the Coulomb potential on 
the wave functions is neglected and neutron and proton wave functions are assumed to be 
the same. 



B. Choice of two-body interactions. 

We have first used a zero range density dependent interaction. The general form of such 
a force writes as: 

V'(r..r,) = -V„ {l - „ (f<^^l±l2mYh(r, - r,) (14) 



With the parameters of Garrido et al. HI], Vo=500 MeV.fm^, q;=0.47 ,?7=-0.1 we get 



much too weak binding. These parameters were fitted so that to reproduce the gap in 
nuclear matter calculated using the Paris potential, then they do not take account of the 
medium effects which are expected to be important. Therefore we have proceeded in a 
different way. We first use a density independent zero range neutron-proton interaction 
where the strength Vq is fitted to give the binding energy of deuteron. The calculation leads 
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to a relation between the T=0 and T=l strengths given by 

Vo{T = 0) = VoiT = 1) |l - (^::pLy\rctg (^)'^'} (15) 

where is the deuteron binding energy and ec the cut-off on the nucleon energies which in 
a nucleus should be counted relative to the bottom of the average one-nucleon potential. In 
our calculations we take for a nucleon in the field of our ^^C-core a cut-off of 10 MeV which 
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is equivalent to ec=60 Mev and corresponds to: 



Vo{T = 0) = 1.36Vo(T = 1) 
= 682MeV.fm^ 



(16) 
(17) 



A similar ratio between the strengths of the T=0 and T=l pairing interactions has been 

n 

found by Satula and Wyss [13|] in their study of even-even nuclei. However a zero range force 
can be considered only as a simple substitute to a more realistic one, and before introducing 
medium effects in our calculations we have tested this substitution and made the calculation 
with the Minnesota free nucleon-nucleon force This force has one short range repulsive 
components and two long range attractive ones with all exchange terms and has been fitted 
on nucleon-nucleon scattering lengths and deuteron binding energy. The results with the 
two forces show very similar results with the same underestimation of two-body correlations 
suggesting the necessity to take account of medium effects on the neutron-proton interaction. 
Indeed when two nucleons are added to a core their mutual interaction will be modified by 
the presence of the other nucleons. Strong two-body correlations which manifest themselves 
as very collective low energy vibrationnal states can induce a modification of the interaction 
through an exchange of these collective phonons between the pair of nucleons. Such phonon 
exchange contribution to the two-neutron pairing force have been studied by Barranco et al. 
[l^ and found to be responsible for about half of the gap in the isotopes ^Ca, "^Ti and "^Sn. 
In our case the ^^C core has a low 2+, T=0 state at 4.4 MeV with a very strong collective 
transition amplitude /?2=0.6 and a less collective 3^, T=0 state at 9.6 MeV with (3^=0.4 
[isj ]. We can expect that these two states will give most of the effect and we include both 
of them in our calculation. 

The diagrams corresponding to the exchange of phonons are represented in Fig.l for the 
three types of matrix elements entering in the RPA equations. The diagrams a) -b) concern 
the matrix elements < a|V|6 > while the diagrams c) -d) and e) -f) concern respectively 
the matrix elements < a|l^|a > and < > of eqs.(5-6) . The calculation of their 

contribution to the RPA matrix elements is given in the Appendix when each vertex of the 
diagrams is replaced by a phenomenological anzats. These matrix elements which have to be 
added to the bare matrix elements depend on the eigenvalues of the RPA equations which 
therefore are now nonlinear equations and will be solved by iteration. 

Once one has seen the effect of phonon exchanges on the ^^N and spectra, the last 
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a) 



b) 



d) 



e) 



f) 



FIG. 1: Diagrams corresponding to a one phonon exchange between neutron-proton pairs appearing 
in the RPA equations. 

step is to look if a phenomenological effective interaction such as the Gogny interaction may 
be interpreted as a bare interaction corrected by the diagrams of Fig.l. In the T=0 channel 
the Gogny interactions, Dl [2] and DIS ^, have a density dependent component of zero 
range. As any zero range interaction, it diverges and its fitted strength depends on the 
cut-off of the single nucleon energies. Since our cut-off for neutrons and protons is quite low 
(10 MeV) we have used a slightly increased strength, as will be discussed below. 



III. RESULTS 

Our configuration subspace is restricted to neutron states up to 10 MeV and to the 
corresponding proton states. 



A. 0+, T=l states 

Previous calculations jl| have shown that the 0^, T=l states in the Li, Be and C isotopes 
are very well described in the two-neutron RPA model with the Gogny forces or their zero 
range density dependent substitute ^. ^^N and ^''B are the analogs of ^^C-^°C where the 
two- neutron pair is replaced by a neutron-proton pair added or substracted from a ^^C- 
core. Because of charge invariance of strong interactions the same T=l effective interactions 
should be able to describe the two kinds of nuclear systems. The calculation has been 
performed for ^^N and ^''B and yields to the results reported in Table II for the zero range 
density dependent force and the Gogny interactions Dl and DIS. We see that the three 
forces give close results as already found in ref[l] for ^^C and a good agreement with the 
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a) b) c) exp.) 




-10.12 -10.18 -10.03 -10.14 
-3.28 -3.10 -3.06 -3.8 


10b 


29.7 29.0 29.24 29.15 



TABLE II: Energies(in MeV) of the 0~*',T=1 states with respect to the ground state of ^^C in ^^N 
and ^'^B obtained with: a) the zero range density dependent force, b) and c) the Dl and DIS 
Gogny interaction respectively. In the last column are given the experimental energies. 
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181 ] even though the excited 0^ state is slightly too high as it 



measured energies 

was in ^^C. We may conclude that we have a good test of the validity of our two-body RPA 
model as well as a further confirmation of the efficiency of the Gogny effective interactions 
to describe light nuclei for which it was not designed. 



B. 1+ and 3+, T=0 states 

All results of this section are presented in Tables III-V and Figures 2 and 3 for ^^N and 
. 

With the zero range interaction of eq.(14) and the parameters of ref.[ll], Vo=500 
MeV.fm^, 1] =-0.1 and a=0.2, our results desagree with experimental spectra. In partic- 
ular the lowest 1"*" state in ^^N is above the lowest 0"*" T=l which is then the ground state 
in desagreement with experiments 16 1. 

Instead we have first made the calculations using the bare density independent zero 
range neutron-proton interaction with the strength Vo=682 MeV.fm^ of eq.(17) and with 
the Minesota force [sl. The energies of the lowest states referred to the theoretical ground 
state energy of ^^C as defined in eqs.(7-8) are given in the table III. We see that the energies 
are very similar for the two bare interactions but that there is a significant desagreement 
with experimental energies for both nuclei. The desagreement is still more pronounced for 
where the ground state is found as a 1+ state while experimentally it is a 3+ state. 
As expected we see clearly that a free neutron-proton interaction is not sufficient to give 
the binding of two nucleons inside a nucleus. Then we have introduced medium effects 
due to the exchange of phonons as explained in section II. B and in the Appendix and have 
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a) b) c) d) e) exp. 


14n 1+ 
1+ 

3+ 


-11.7 -11.9 -12.4 -12.8 -12.5 
-5.8 -4.61 -5.75 -6.9 -6.9 -6.3 
-3.26 -4.27 -5.55 -7.25 -7.2 -6.05 


3+ 


27.6 28.7 28.1 29.1 28.4 28.1 
30.3 29. 28.6 28.3 27.7 27.4 



TABLE III: Energies of the lowest 1^ and 3^,T=0 states in N and B with respect to the ground 
state energy of ^^C obtained with different interactions: a) the bare zero-range interaction -b) the 
bare Minnesota interaction -c)the Minnesota force plus the exchange of phonons -d) and e) the 
Gogny force with respectively the experimental and calculated IP3/2 proton energy . In the last 
column are given the experimental energies. 

calculated the induced matrix elements according to eqs. (21-29). The results are again very 
similar for the two interactions and we show in the tables and figures those obtained with 
the Minnesota interaction. The excitation spectra are shown in Fig. 2 for ^^N and Fig. 3 for 
for the bare and bare plus induced interaction. At the bottom of the figures are given 
the separation energies of a neutron-proton pair in ^^N and ^^C which are directly related 
to the lowest RPA energies in ^^N and respectively. 

We see an improvement for the ^^N excited spectrum as well as for the neutron-proton 
separation energies. The excitation spectrum of ^^N (see Figure 2) shows now a quite 
good agreement with the experimental spectrum. Since we take in the definition of our 
neutron-proton subspace nucleon states up to 10 MeV we have in the higher part of the ^^N 
spectrum a large number of states but we have given only those called pairing vibrational 
states which have amplitudes on several neutron-proton configurations. The other levels 
have a smaller probability to be seen experimentally |l9|. For these high energy levels we 
have not calculated the contribution of phonon exchange. Our code is inefficient when we 
have several very close eigenstates but we expect the induced matrix elements to be weak 
because of energy denominators. The contribution of phonon exchange is very important for 
the low energy states and improves significantly the energy spectrum where now the ground 
state energy, or equivalently the neutron-proton separation energy, is very good. At about 
6 MeV we can reproduce the experimental 1+- 3+ doublet with the right order of levels 
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„,, 2 levels ^ 2 levels ^ 3+ 



L]^+ 

S„pC^''N)=ll,9 McV 12,4 MeV 12,8 iVlcV 12,5 Ml^V 

»} b) c) 

FIG. 2: Excited spectrum of ^'^N calculated with: a) the bare Minnesota interaction-b) the same 
with the exchange of phonons-c) the Gogny interaction ( DIS with t3=1600 MeV.fm^. The exper- 
imental spectrum is given in the last column. Below are given the ccalculated and experimental 
neutron-proton separation energies. 



- - 1+ 

L 1+ „ " . 3+ 



I 

LI 



S^^(^2 C)=28,7 MeV 28,1 MeV 28,3 MeV 27,4 MeV 

a) ^ c) d) 

FIG. 3: Same legend as Fig. 2 for ^'^B . The numbers at the bottom of the figure are the neutron- 
proton separation energies, calculated and experimental, in ^^C. 

and energies very close to the experimental ones. The main RPA amplitudes are given in 
Table IV for the lowest states . We see that the RPA Y-amplitudes with two nucleons in 
the lp3/2 shell are not negligible what indicates the presence in the wave function of the 
^^C ground state of configurations with at least two holes in this shell, in agreement with 
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E*(MeV) 


(lPl/2)' 


(1^5/2)' 


(2s)2 


(2s,ld5/2) 


(lp3/2)2 


1+ 





0.88 


0.24 


0.1 




-0.24 


1+ 


6.6 


-0.22 


0.50 


0.74 




-0.11 


3+ 


6.8 




0.50 




0.78 


-0.13 



TABLE IV: RPA amplitudes in ^^N 





(lPl/2)?+,0 


(lP3/2)?4 


-,0 (lP3/2)3+,o 


Min 


-3.42 


-3.74 


-5.19 


Mill' 


-3.42 


-4.03 


-5.49 


PWBT 


-3.45 


-4.16 


-6.04 



TABLE V: Diagonal matrix elements (in MeV) for the two nucleons in the lpi/2 or shells 
coupled to J^=l+ or 3+ calculated with the bare Minnesota interaction (Min), when the phonon 
exchange contribution is added (Min') and those fitted in ref.[24] by Warbuton and Brown(PWBT). 



the shell model calculation of Cohen and Kurath for example [20| . In the higher part of 
the spectrum, above 10 MeV, we get a group of two 1+ and two 3^ states which may be 
identified with the same experimental group. However we miss the experimental 1"*" level at 
3.95 MeV which is very likely formed mainly of a ^^C core excited to its 2+ state at 4.4 MeV 



coupled to the n-p pair in its 1+ ground state. 



2l| and ^^0(7,np)^^N reactions 
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his is also suggested by the analysis of C 



231]. RPA is not able to describe such state 



(6Li,a)i4N 

since it relies on the assumption of an inert core in its ground state. In the same way we 
miss the 1^ and 3^ states at about 9 MeV which again are very likely due to the coupling 
of the same 2^ in ^^C with the neutron-proton pair in the 1~^ or the 3"*" excited states at 
about 6 MeV. This assumption is supported by the presence in the experimental spectrum 
Qof two T=0 levels, a 2+ at 8.98 MeV and a 5+ at 8.96 MeV which very likely belong to 
the same multiplet. 

In the results are improved by phonon exchanges but not enough to get the 3"'" state 
as the ground state. The energy of the 1^ when referred to the core ground state energy, is 
close to the experimental value but the 3"*" is too high and appears again as an excited state. 
This 3+ state is a nearly pure (p3/2)^^ state weakly affected by RPA correlations because 
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the lowest possible configuration with unoccupied neutron and proton unperturbed states 
(anomalous RPA configuration) is the (2s,ld5/2) configuration, much higher in energy. The 
correlated energy is therefore mostly due to the diagonal matrix element ((^3/2)^1^1(^3/2)^) 
which appears to be too weak as is discussed now. Indeed we compare our matrix elements 
with those of Warburton and Brown [2JI which are determined by least-squares fits to 216 
levels in A=10-22 nuclei. They use a (lp,2s,ld) shell model space, compatible with our 
subspace, at least for the lowest states. In Table V are reported their fitted matrix elements 
and ours calculated with the Minnesota plus phonon exchange when the n-p pair is in the 
Ip-shell. We show the diagonal matrix elements for the n-p pair in the lp3/2 and lpi/2 shells 
coupled to spin I"*" and 3''' which are for a large part responsible for the binding energy of 
the lowest 1"^ states in ^^N and and the 3"'" state in . The same comparison has been 
made with the matrix elements fitted by Cohen and Kurath {2^ in a smaller model space but 
the conclusion is qualitatively the same. We see a good agreement when the two nucleons 
are coupled to J=l^ but a too weak binding for J=3~''. It shows that something is missing 
in our force and a possible reason why it is too weak could be found in more recent shell 
model calculations. Indeed in a no-core shell model and with the Argonne V8' two-nucleon 
force , Aroua et al [26] get a good description of the 1^ ground state of ^^N while Caurier 
et al.[25| and Navratil and Ormand [9] show that with the same force and the same model 
it is not possible to get the ground state spin of ^°B . The authors of ref.[9] show that it 
is correctly obtained only when they include the Tucson-Melbourne three-body interaction. 
The necessity to use a three-body force in the description of ^°B is also clear from quantum 



Monte Carlo shell model calculations of B 
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A genuine three-body force is out of the 



scoop of the RPA since we work with two nucleons only. However this three-body term could 
induce an effective two-body contribution which is not included when we use the Minnesota 
force plus phonon exchange. Intuitively we expect this term to be more important for two 
lp3/2-nucleons inside the ^^C, therefore for the description of ^°B , where they will interact 
with a third nucleon off the six other nucleons in the same shell than for the two nucleons 
added in higher shells to describe ^^N which will have less interaction with a third nucleon 
inside ^^C. This might explain why our calculations failed in ^°B but are very satisfying in 

At last we have made the same RPA calculations with the Dl and DIS Gogny forces 

aa. 

The two forces give the same results and we discuss only those obtained with DIS. In the T=0 
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channel this force has a zero range density dependent component with a strength t3=1390.6 
MeV.fm^. Any zero range interaction is divergent and the fitted strength depends on the cut- 
off on single-particle energies. Since our cut-off is somewhat low we have increased t3 in order 
to get close to the measured energy of the 1^ ground state in -"^^N . With t3=1600 MeV.fm^, 
a value slightly stronger than the genuine value, we get a binding energy for the n-p pair of 
-12.7 MeV instead of -12.5 MeV, the experimental value. In the third columns of Figs. 2 and 
3 we show the excited spectra obtained for ^^N and respectively. We see that for ^^N 
the levels are in very close agreement with experiment and very close to those obtained with 
the Minnesota interaction plus the phonon exchange contribution with however an inversion 
of the 1+- 3+ states at about 6 MeV. This result suggests that the Gogny interaction, 
often thought as a G-matrix, includes implicitly the phonon exchange contribution. In 
the 3"*" state is obtained as the ground state as it should, contrarly to what we got 
previously with an effective two-body interaction. The 1"*" state is now the first excited state 
with an excitation energy in agreement with the experimental value and with shell model 
calculations of ref[9]. Therefore according to our previous discussion on shell model results 
we may conclude that, through its density dependent term, the Gogny interaction includes 
implicitly an effective two-body contribution coming from three-body forces. However the 
neutron-proton separation energy in ^^C given as the difference between the ground state 
energies in and ^^C is too high by 0.8 MeV but a still larger overestimation is observed in 
the shell model calculations (see their tables V and VII) . Note that if we use the calculated 
energy for a lp3/2 proton (see Table 1) instead of the experimental energy this discrepancy 
is strongly reduced. 

IV. CONCLUSIONS 

In the framework of a two-particle RPA model applied to the description of and ^^N 
formed of a correlated core of ^^C in its ground state minus or plus a neutron-proton pair, we 
have performed a detailed analysis of the T=0 effective nucleon-nucleon interaction. First we 
have shown that for the 0+ T=l states the results are as good as they were for ^^C-^°C (where 

the neutron-proton pair is replaced by a two- neutron pair) when the Gogny interactions, Dl 
and DIS, or their zero range equivalent are used. The same two-nucleon RPA model has 
been applied to 1"*" and 3"*", T=0 states and an attempt to analyse the contains of an effective 
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interaction has been made. We have shown that an effective interaction constructed as a 
bare interaction, from Minnesota or of zero range, fitted to the deuteron binding energy, 
complemented by medium effects due to phonon exchanges between the neutron-proton pair 
gives a very good representation of the -"^^N levels but fails to reproduce the 3^ ground state 
of . By comparing with shell model calculations wc arc able to suggest that this is due 
to the presence of a three-body component in the interaction which are not included in 
our derivation. This additionnal component will not spoil our good results in ^"^N but will 
improve those in ^'^B . At last our study suggests that the T=0 Gogny interaction which yields 
good agreement with measurements in both nuclei, includes empirically both the effect of 
phonons exchange in the effective interaction and the effect of a two-body component coming 
from the presence of a three-body interaction and included in the density dependent term. 

I am very grateful to E. Viggezzi and F. Barranco for pointing to me a crucial problem 
in the early stage of this work. Many thanks are also due to N. Auerbach and B. Barrett 
for their interest in this work and for a careful reading of the manuscript. 

APPENDIX A: MATRIX ELEMENTS OF THE INTERACTION INDUCED BY 
PHONON EXCHANGE. 

The matrix elements of the induced interaction of Figure I are calculated between two 
neutron-proton pairs coupled to spin J,M and isospin T=0. 

Each vertex in the diagrams is replaced by a phenomenological expression. Let's call 
L,Mi the angular momentum of the phonon, the transition density from a zero-phonon 
state to the one-phonon state is written as: 

< lph\V\Oph UlRo^^Y^'^^u) (A1) 
L dr 

where U(r) is the average one-body potential assumed to be a Saxon-potential such as: 

where the function f(r) and the values of the parameters are given in section 11. The 
collective amplitudes, Pl, can be fitted on the experimental values of the B(EL) or on 
proton inelastic scattering cross sections. 

The two-nucleon wave function is constructed by coupling the two-nucleon states to a 
total spin (JM) as: 

15 



\jij2, JM >^ J2 < 3i32mim2\JM > \lijimi > \l2j2m2 > (A3) 

mim2 

where the nucleons 1 and 2 are either in occupied states, ai 0:2, or unoccupied states, Oi 02 
with energies , or respectively. The total wave function has to be antisymmetric 

so that for T=0 states the spin-space wave function has to be symmetric. Therefore the 
general expression for the antisymmetrised matrix element of the induced interaction, V^^^, 
is obtained as: 

(A4) 



L 



with: 



1 /?-r^- I 



jl J3 ^ j ( ^2 J4 ^ I I Jl J2 J 

0.5 -0.5 J I 0.5 -0.5 ) I j4 Ja ^ 



> < l\g\3 >< 2\g\4: > (A5) 



VeL = 
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J1J2J3J4 
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i2 is 
0.5 -0.5 
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32 


31 


> 
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II 






> 




U 


33 


L 

> 



< m\g\n > 



0.5 -0.5 

L 9{rWir.3mi^)<t^injn{r)r^dr (A7) 

J u 

where the 0's are the one-nucleon radial wave functions. 

While the quantities V^l and V^l have the same expressions whatever are the initial and 
final two-nucleon states, occupied or unoccupied, the expressions of D^^ and D^.^ which 
come from energy denominators have to be derived in three different cases corresponding to 
the diagrams of Fig.I: 

If (12)=(ai a2), (34)=(a3 a4) (diagrams a) and b)) 

DdL = (e„, + e„3 + ujl)Y^ + [Q - (e„, + e„, + u^l)]'^ 
If (12) = (ai a2), (34) = (q;i 02) or the inverse (diagrams c) and d)): 

DdL = -[Caa - + " " ^ai + ^^l]"^ 

DeL = -[Cai - + ^l]~^ - ^2 " ^ai + ^^l]"^ (AlO) 



(A8) 
(A9) 
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It is straightforward to show that: 

< aia2|l^md|QiiQi2 >=< Oiia2\Vind\aia2 > (All) 
If (12)=(ai 02), (34)=(q;3 a^) (diagrams e) and f)): 

DdL = -[^ - Caa - fas + - - - 6^^ + UJl]~^ 

DeL = -[O - - + (jJl]~^ -[n- e„2 - e^^ + 0;^,]"^ (A12) 

In these equations Q is the eigenvahic of the RPA eqs.(5-6) and uj^ is the energy of the 
phonon L. When we add these phonon exchange contributions to the bare matrix elements, 
the RPA equations become non-linear and will be solved by iteration. 

Note that the approximation of eq.(18) implies that one can defined an equivalent two- 
body term which has to be added to the bare neutron-proton interaction which is the sum 
of separable terms of the following form: 

SV = ^^pP-Y^{u;)Y^{u;') (A13) 
dr dr' ^ ^ ' ^ ^ ' ^ ' 
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